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We study the diffusion-limited reaction A + A «-» A in various spatial dimensions to observe the 
effect of internal fluctuations on the interface between stable and unstable phases. We find that, 
similar to what has been observed in d = f dimensions, internal fluctuations modify the mean-field 
behavior predictions for this process, which is given by Fisher's reaction-diffusion equation. In d > f 
the front displays local fluctuations perpendicular to the direction of motion which, with a proper 
definition of the interface, can be fully described within the Kardar-Parisi-Zhang (KPZ) universality 
class. This clarifies the apparent discrepancies with KPZ predictions reported recently. 



PACS numbers: 5.40.-a, 64.60.Ht, 68.35.Ct 

Fluctuations in the macroscopic behavior of reaction- 
diffusion (RD) systems could play an important role 
whether these fluctuations are produced by the discrete 
nature of the elementary constituents (internal) or by 
environmental random variations (external) . Often both 
type of fluctuations are neglected in the theoretical treat- 
ment, thus describing the system by an action-mass type 
equation (mean- field approximation). Fluctuations in 
RD systems can, for instance, give rise to instabilities |Q , 
modify the reaction front velocity [|]]|] i a U°w the system 
to reach new states absent in the mean-field description 
, or produce spatial correlations in the system which in 
turn can dominate the macroscopic system behavior ||. 
While these effects appear in different situations, there is 
a theoretical and experimental interest in the problem of 
front propagation in RD systems. The simplest example 
is that of invasion of an unstable phase by a stable one, 
which at the mean-field approximation level is described 
by Fisher's equation 

^ = DV 2 p + k lP -k 2 p 2 (I) 
at 

where p(x,t) is the local concentration (x 6 R d ) char- 
acterizing the system. Equation ([!]) arises in the macro- 
scopic description of many processes in physics, chem- 
istry and biology and is a generic model for reaction front 
propagation in systems undergoing a transition from a 
marginally unstable (p = 0) to a stable (p oq = fei/fe) 
state. Thus, for initially segregated conditions, i.e. 
p(x) = p cq for xu < and p(x) = for Xu > (with 
x = (xx,Xn), x± £ M. d± , X|| € K and d = d± + 1), the 
solution of (Q) is a front invading the unstable phase and 
propagating along Xu with a constant velocity v > v m i n — 
2(kiD) 1 / 2 which is selected depending on the initial con- 
dition according to the "marginal stability criterion" . 
At the same time, the front broadens until it reaches a 
finite width f ~ 8(L>/fci) 1/2 - 

The question of how faithfully continuum equation (Q) 
resembles the macroscopic front dynamics of microscopic 
RD discrete systems has drawn a lot of attention recently. 
In particular, much attention has been devoted to micro- 
scopic stochastic models like A + A A or A + B ^ 2 A, 



where A and B are active species. Discreteness of those 
systems is responsible for fluctuations in p(x,t) and in- 
troduces an effective cutoff in the reaction-rates which 
modifies the properties of the front. Most of the studies 
have concentrated on observing how the microscopic sys- 
tem approaches the macroscopic behavior described by 
Eq. (Q) in d = 1 when N — > oo, where N is the num- 
ber of particles per site Q . Using van Kampen's system 
size expansion || or field-theory techniques S , the meso- 
scopic dynamics of the microscopic Master equation can 
be expressed in terms of a Langevin equation which, in 
the case of the A + A <-> A scheme, reads [jl(J 

dp 



dt 



DV 2 p + kip - k 2P 2 + y/(k lP - k 2 p 2 ) V (x, t) (2) 



where n(x, t) is a uncorrelated white noise, {77(05, t)) = 0, 
(?7(£c,<)?7(a; / ,t')> = 2N~ 1/2 6{x - x')S{t - t'), and N is 
a reference level of number of particles, which is pro- 
portional to N. In the limit N — * 00, the noise term 
in (0) (which reflects the fluctuations in the number of 
particles) vanishes and we recover ([!]). In this limit, the 
effective cutoff in the reaction term imposed by the dis- 
creteness seems to be the leading contribution, and it is 
possible to derive corrections to the velocity which yield 
the well known result v — v m i n ~ 1/ log 2 N [p 10 1. 

The purpose of this paper is to extend this analy- 
sis of discrete microscopic models to higher dimensions. 
In d > 2 the front position is given by an interface at 
xii = h(x±,t) which separates the stable from the un- 
stable domain. Due to the microscopic fluctuations one 
expects the interface to roughen and its fluctuations to 
be described asymptotically by the Kardar-Parisi-Zhang 
(KPZ) Jll[] equation, which features the simplest and 
most relevant non-linearity (in the renormalization group 
sense) that considers local and lateral growth: 



^=v + DVlh+±(V ± hy< 



+ V2c^ri ± {x ± ,t), (3) 



where r)j_(xj_,t) is an uncorrelated white noise, and Vi 
is the divergence operator defined over the substrate x± . 
The interface is described by its mean position h(t) and 
its roughness, 
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w 2 (L ± ,t) = ([h(x ± ,t)-h(t)} 2 ), (4) 

where (• • •) means average over different realizations and 
the bar denotes average over the substrate x± € 
(L± is substrate lateral length). As it is well known, in 
the KPZ equation w obeys a scaling form w(L±,i) — 
tPf(t/L z ± ), where w ~ t p if t < L z and w = w sat - L Q 
if t 3> L z , where a = flz is the roughness exponent [fL2[. 

However, studies of microscopic realizations of Eq. (Ill) 
have questioned the applicability of (9) to describe the 
interface fluctuations in front dynamics. Riordan et al. 
|Q studied the A + A «-> A scheme in various dimen- 
sions concluding that the interface roughens in time for 
dx < dx,c = 2 with /3 d±=1 = 0.27 ± 0.01 and f3 d±=2 0. 
These values strongly differ from those of KPZ, which 
are f3 d±=1 = 1/3 and d± = 2 = 0.245 @. For d± > d x , c , 
they observed that the interface is flat and described 
by Eq. (|l|). However, the determination of these val- 
ues is hampered by the method used to obtain them, 
namely, using the front width of the projected density 
p(x|i , t) obtained by integrating out all perpendicular 
dimensions x^_. To overcome this problem, Tripathy 
and van Saarloos JUm studied the interface dynamics in 
d± = 1 in terms of the coarse-grained density of par- 
ticles (see below) and obtained Pd±=i = 0.29 ± 0.01, 
a d± =i = 0.4 ± 0.02. The authors interpreted this de- 
viation from the KPZ exponents (otdj_=i — 0.5) as an ef- 
fect of the parallel dimension, which cannot be integrated 
out to obtain an effective interface description in dj_ di- 
mensions and thus, the exponents correspond to those of 
the KPZ equation in d± + 1 dimensions (a d± =2 — 0.393) 
, although the interface is a surface defined in a d±_- 
dimensional substrate. More recently fllq] , it has been 
conjectured that this situation is more general and oc- 
curs in the dynamics of the so called "pulled fronts" in 
which the region ahead of the front is important to de- 
termine not only its deterministic properties (like the ve- 
locity) but also the fluctuations of the effective interface 
dynamics |Tq] . Thus, it is important to check the validity 
of this conjecture for d± > 1. Although our results are 
consistent with those in p5| , a closer analysis indicates 
that there is a crossover to the asymptotic behavior of 
KPZ equation in d± dimensions. 

Our simulations are done on a x Ln lattice with 
Lii 3> L±. At each point of the lattice we define the 
occupancy number n(x,t) to be 0, 1 depending whether 
the site is occupied by a particle or not. At each time 
step a particle is chosen and can perform any of the pro- 
cesses in Fig. |l|a which, in the mean-field approximation, 
yield k\ = p and k 2 = 1 + p, i.e. p c(l — p/(l + p) 
ipHll . Diffusion rate is chosen to be D = l/(2d). Ini- 
tial conditions are set up by distributing particles uni- 
formly over a domain x with density p cq and 

(i) 

i|l of the order of hundreds of lattice spacings. The 
interface is determined using a local coarse-grained den- 
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FIG. 1. (a) Microscopic rules for our model: diffusion 
(left) and reaction (right). The model is isotropic, i.e., all 
movements are equally probable in any spatial direction, (b) 
and (c) Snapshots of the stable domain for fj, = 0.1 (b), 
fj, = 0.5 (c) and I — 5. 

sity pi(x,t) = l~ d J2seD l n ( s 't) wnere Di is a squared 
domain of lateral size I around the x site. The stable 
domain of the front is defined by those sites for which 
pi(x,t) > 5Pcq and the position of the interface h(x±,t) 
is the largest value of asn inside the stable domain for 
fixed x± g5). 

Simulations show that the interface advances linearly 
in time, h = vt, where from we can extract the velocity of 
the front in any dimension. Values obtained are reported 
in Fig. |^ and compared with the deterministic prediction 
v = Umin = Is/Dp, of Eq. ([!]). We recall that in d = 1 
the stochastic model admits an exact solution [ fl7|| , and 
the front advances with velocity v — Dp. The differ- 
ence between both expressions reflects the importance of 
fluctuations when N = 1 and d = 1, and agrees with 
the breakdown of the approximation v — u m i n ~ log -2 N 
when N — 1. Interestingly, for d > 1 we recover the 
law v ~ rt 1 / 2 , although the value is still corrected due to 
the fluctuations in the number of particles. Nevertheless, 
the velocity change is smaller when the dimension is in- 
creased, as fluctuations get averaged over the transverse 
dimensions. Moreover, assuming that Eq. (§) holds, we 
can obtain the value of A by tilting the substrate with an 
overall slope m and measuring the velocity. As expected 
we obtain v(m) = v(0) + \v{Q)m 2 , i.e. A = v(0) @. 

For d — 2 the interface roughens subject to internal 
fluctuations, and the roughness grows like w(t) <~ t 13 
until it saturates to a constant value u> sa t ~ (see 
Fig. |). When p = 0.1 we obtain (3 d±=1 = 0.27 ± 0.01 
and a — 0.40 ±0.01, consistent with Neverthe- 
less, a closer inspection of the data suggests that there 
is a crossover from non-Gaussian fluctuations at small 
scales to Gaussian ones at large scales. Although this 
crossover could be guessed, for instance, from the value 
of w sat (L±) for p — 0.1 (which is seen to change its 
slope in a log-log plot), it becomes clearer if we look 



2 



10 ' 



10" 




10' 



10° 



10' 



10 



10" 



FIG. 2. Front velocity for different dimensions. Dashed 
lines are mean-field prediction v = 2(D/i) 1//2 and d = 1 exact 
solution v = Dp. Values for d > 2 are divided by 2D 1//2 . Er- 
ror bars are smaller than the symbol size. Inset: renormalized 
velocity as a function of the tilt m for d± — 1. Dashed line is 
the parabola m /2 and wo = v(m = 0). 



at higher moments of the h(x±,t) distribution, e.g. the 
skewness, S = ((h(xj_,t) — h) 3 ') /w 3 ^ 2 (t) and the kurto- 

sis, K. = ((h(x± } t) — hY)/w 2 (t) — 3 when t — > oo as 
seen in Fig. |[ Both approach the Gaussian asymptotic 
regime given by Eq. (|^) (i.e. S = K, = 0) when L± — > oo. 
Non-Gaussian fluctuations in the system are due to the 
vague definition of the interface for small values of p and 
I. In Fig. [l]b,c we show different snapshots of the sta- 
ble domain calculated through the coarse-grained den- 
sity, pi{x). It is apparent that when the number of par- 
ticles inside the domain Di is small, the definition of the 
interface is vague and implies large steps (or large slopes) 
which give non-Gaussian fluctuations in h(x,t) at small 
length scales. Nevertheless these intrinsic fluctuations of 
the interface are bounded and restricted to small length 
scales. The measured roughness is naturally decomposed 
into contributions due to intrinsic fluctuations, Wi, and 
long- wavelength fluctuations wq fl8|,|l9|| 



w 2 (L,t)=w 2 +w 2 .(L,t), 



(5) 



where u>i is a non-universal constant that depends on 
the number of particles, n;, inside the domain Di used 
to determine the interface. Direct fit of expression (|^) 
to a power law when t — > oo gives lower values of the 
exponent a, while we recover a = 1/2 when fitting the 
data in Fig. |^ using expression (|^) and assuming that 
wo ~ L a , showing that long- wavelength fluctuations are 
indeed described by the KPZ equation in d± dimensions. 
In order to reduce the effect of the intrinsic width, we 
can enlarge the size of the domain and/or the value of 
/x, as we can see in Fig. [|, where a direct fit to a power 
law of the data for p = 0.1 and I = 9 gives a = 0.47 ± 
0.01. In this sense, ni plays the role of a noise-reduction 
parameter. The situation then is completely reminiscent 
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FIG. 3. (a) Roughness as a function of L± for d — 2. 
Dashed lines correspond to power laws with a = 0.4 and 
a — 1/2. (b) Same as in (a) but for the directed percola- 
tion model with c = 0.1. (c) Skewness (circles) and kurtosis 
(squares) as a function of L± for /i = 0.1 (full symbols) and 
p = 0.5 (open symbols). Error bars are about 0.05 but have 
been omitted for clarity reasons. 

of what happens in the Eden model [pi and has been 
also observed in many other RD systems jl9| . 

In higher dimensions results are analogous. In d = 3 
we do observe a significant increase in the roughness with 
L±, contrary to what was observed in |p"3f . The reason 
for this discrepancy is that the roughness measured in 
|]l3| is also a linear combination of the long-wavelength 
roughness wo(L±,t) and the front width, £. In d = 2 
we have wo(L±,t) > £, while in d = 3, Wo(Lx,t) ~ £ 
(for the simulated values of L±), and this is responsi- 
ble for an apparent value of ad ± =2 — 0. Actually, us- 
ing the structure factor S(q±,t) — (h(q±,t)h(—q±,t)) 
where h(q±,t) is the Fourier transform of h(x±,t), we 
observe in Fig. a crossover to the asymptotic scaling 



S(q± « l,i -> oo) 



with otd ± =2 — 0.393 of 



KPZ equation [oscillations for q > 1 arc due to artificial 
correlations in h(x,t) introduced by the coarse-grained 
density] . Finally, in d = 4 we observe two different behav- 
iors: for small values of // the roughness seems to saturate 
to a value independent of L± while for = 0.5 we recover 
the asymptotic prediction of Eq. (||), ad ± =3 = 0.313 [ pT| . 
This difference may be due to the existence of a phase 
transition for the KPZ equation in dj_ > 2, in which, 
depending on the control parameter g 2 = \ 2 aD~ 3 , the 
interface is flat (for g < <? c ~ 6.28) or rough (for g > g c ). 
Our results suggest that the interface is flat for p = 0.1 
(i.e. w 2 does not scale with L±) while it is rough for 
fj, = 0.5, although more simulations are needed. 

We have also performed simulations for the N = 1 di- 
rected percolation model in which the decay A — > is 
allowed with rate c. The model presents a non-trivial 
phase transition for finite values of fi and c in which the 
stable state changes from an active (p eq ^ 0) to an ab- 
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'kip rj(x,t), 



(6) 



FIG. 4. Structure factor as a function of the wave number 
q for d — 3 and Lj_ = 128. Dashed line correspond to the 
KPZ prediction in d ± = 2, S(q) ~ q~ 2+2ct with a = 0.393. 
Inset: Saturation roughness as a function of L± in d = 4. 
Dashed line corresponds to the power law L°f with a = 0.313. 

sorbing state (p cq = 0) for d < 4 j20|. In this case, the 
mesoscopic dynamics of the system are given by [H,^| 

dp 
dt 

the only difference with @ being the form of the noise 
term. Away from the transition point (which occurs 
when p — 0.22 ±0.01 for fixed value of c = 0.1), interface 
fluctuations behave similarly to the previous model (see 
Fig. |j^b), confirming that our results do not depend on 
the specific choice of the microscopic rules, as long as the 
mean- field limit of them coincides with (|l]) . 

In summary we have studied two microscopic stochas- 
tic realizations of Fisher's equation with TV = 1. It has 
been shown that with a proper definition of the interface, 
long-wavelength fluctuations are fully described within 
the KPZ scaling in dj_ substrate dimensions. Thus, for 
d < 4 there is not an upper critical dimension above 
which (0) is valid. Only for d — 4 does the interface 
seem to undergo a phase transition between a flat inter- 
face [which can now be explained by (Q)] and a rough 
phase, depending on the microscopic parameters. 

Finally, we comment on the possibility to observe the 
conjecture in fl6| l using microscopic RD systems. It is 
possible to argue that the results in this paper deviate 
from those of the conjecture because the limit N — > oo 
in which Eq. (|l|) is valid has not been reached. Thus for 
large values of N < oo |H| we should be able to observe a 
crossover to the universality class proposed in [|l6| . How- 
ever, the conjecture requires that the noise term scales as 
pr\{x, t) when p — > 0, which we call multiplicative noise. 
This type of noise can be obtained, for instance, assuming 
that the reaction rates fluctuate in time due to environ- 
mental variations (external noise). As argued in ]22] ], it 
is impossible to find such an internal noise in a micro- 
scopic RD model. For example, in the models analyzed 



here, the noise term scales as p 1 ^ 2 rj(x,t) for p — > [see 
Eqs. (|) and (§)]. Thus, even if the parallel dimension is 
relevant to determine the interface properties of the front 
for large values of A, we expect internal fluctuations to 
give a different universality class from the one proposed 
in g6|. 
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Note added: After completion of this paper we became 
aware of a recent pr eprint, R. A. Blythe and M. R. Evans, 
cond-mat/0104395 , in which a similar crossover is found 
for a model in the directed percolation universality class 
of (§) in dx = 1. 
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